Phase signal definition for electromagnetic waves in X-ray crystallography 
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Diffracted X-ray waves in crystals have relative phases regarding the mathematical format used 
to describe them. A forward propagating wave can be denned with either negative or positive time 
evolution, i.e. k ■ r — ut or tut — k ■ r. Physically measurable quantities are invariant with respect 
to the choice of definition. This fact has not been clearly emphasized neither extensively explored 
when deriving well-established equations currently being used in many X-ray diffraction related 
techniques. Here, the most important equations are generalized and consequences of conflicting 
undertaken definitions discussed. 



I. INTRODUCTION 



Electromagnetic waves are in general described by an 
expression such as 



D(r,t) =D (r)e 



iS(u)t-k-r) 



(i) 



where S = ±1 stands for the global phase signal defini- 
tion of monochromatic forward propagating waves. Both 
choices of phase signal (+ or — ) are allowed and both 
can be founded in several theoretical approaches of the 
X-ray diffraction phenomenon in crystals. 

Physically measurable quantities are invariant regard- 
ing two distinct and independent choices: one is the 
choice on how to describe the forward propagating X- 
ray wave, represented here by the options S = — 1 and 
S = +1; and the other is the choice of origin for the 
spatial position r. In the l iterature of X-ray diffraction 
(see for instance Refs.Q|3300) these two dist inct choices 
have been of particular relevance for deriving important 
mathematical expressions, such as those of the structure 
factor of a crystal unit cell and of the atomic form fac- 
tor. Whatever is the chosen option, S — ±1, it can be 
compensated by the choice of origin so that the expres- 
sions are obtained in their standardized format without 
showing any dependence to the phase signal option. 

This standard procedure links the choice of origin to 
the choice of phase signal, which is unreal since both 
choices are independent from each other. The real fact 
is that the available expressions are in agreement with 
only one value of S for each given choice of origin. Con- 
sequently, the propagating waves in the medium must be 
described according to the implicitly undertaken choices, 
otherwise measurable features such as the positioning of 
standing wavefields — formed inside crystals undergoing 
diffraction — as well as invariant phase values — acces- 
sible via 3-beam diffraction experiments^ — will be af- 
fected by conflicting definition of the global phase signal. 

In this work the conflicting points are demonstrated 
and commonly used expressions in X-ray crystallography 
are generalized for any preference of signal choice S. Sig- 
nal shifts owing to geometrical factors are avoided by 
choosing the origin further way from the diffracting vol- 
ume. The consequences only on the choice of S is then 
exploited and discussed in the rest of the paper. More- 




FIG. 1: The relative phase of the scattered wave from R, 
regarding an arbitrary choice of origin at O, is given in terms 
of incident ko and scattered k wavevectors since kMO — — ko- 
r and kON = k ■ r, as demonstrated in Eq. ||5J|. For sake of 
simplicity in analyzing the dependence with the global phase 
signal S, the origin is chosen very far from the charge-density 
distribution p(r), to asssure that (k — ko) ■ r would have the 
same signal, in our case positive, for any point R where p(r) ^ 
0. 



over, to facilitate part of the demonstration, an adapta- 
tion of the Stokes relations for Bragg reflection in perfect 
low-absorbing crystals is presented. 



II. SUSCEPTIBILITIES TO PHASE SIGNAL 
DEFINITION 

A. Structure factor 

As depicted in Fig. 1, the scattered wave 

Afl(fe) = D{B') + D{B) 

towards direction k depends on the relative phase be- 
tween the rays from the origin O and from the electron 
density at any given position r, as for instance at R. In 
the wavefront BB', 

D(B') = AD R e iS[u ' t ~ k{W " +lm7 ' )] 
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and 

D(B) = A£) e^["*- fc (o3+os)] = AD e lv . 

Since the origin has been chosen outside the electron den- 
sity distribution p(r), there is no scattering at the origin, 
i.e. AL>o — 0. Also, the origin was chosen in a such way 
to guarantee that (k — kg) ■ r > to all positions inside 
the volume where p(r) =/= 0, hence 

k(A r R + ~RB') = k(AO + OB) - k(MO + ON) 

and 



coordinates of the atoms in the unit cell, hkl is in fact 
the Miiller index of a given reflection H. 

The ± option in the exponent of Eq. JHJ is not related 
to the choice S but it can be conveniently used to obtain 

Vc H 

Methods for experimental determination of p c (r), or 
structure determination methods, are based on the above 
expression Eq. Q , although it has been standardized for 
the S = +1 choice. 



AD(k) = A DRe iSk(MO+oN) e i<e 

= AD R e lS(k - ko) - r e itp . (2) 

The choice-of-origin phase tp, is common to all elements 
of volume dV at any instant of time, and can be conve- 
niently chosen to provide e i(p = 1. If ADr = D e p(r)dV 
(D e is the scattering amplitude of a single electron), the 
total scattered wave by the volume V is 

D(k) = D e [ p{r)e lS{k - ko ^ r AV. (3) 
Jv 

In X-ray crystallography, the structure factor Fu cor- 
responds to the total scattered wave from an unit cell, 
which is D c (k) — D c Fh /V c (V c is the unit cell volume). 
Therefore, it follows from Eq. © that 

F H = /ne' S2lHT " = \F H \e lSH (4) 

n 

where 2irH = k — ko, H is the diffraction vector of 
reflection H, and /„ is the atomic form factor^ of the 
atom at the position r„ in the unit cell. 

B. Electron density of the unit cell 

The step from Eq. to Eq. (@J is better compre- 
hended when written the electron density of the unit cell 

Pc{r) = y^f n 5{r -r n ) (5) 

c n 

in terms of Dirac <5-functions 

S(r - r„) = e ±j27r - ff - (T " r " ) . (6) 

H 

The sum in H stands for three sum on integer numbers 
h, k and I, running from — oo to +oo since H ■ r n = 
hx n + ky n + lz n and x n , y n and z n are the fractional 



C. Standing waves 

The X-ray reflection coefficient of a crystal with N 
planes (crystal thickness Nd, lattice period d) in sym- 
metric Bragg reflection geometry can be written as 

R N (8) = \R N (e)\e^ H +^ (8) 

where 9 is the rocking curve angle, i.e. the angle between 
the incident X-ray wave and the lattice planes. Sh is the 
phase of the structure factor in Eq. Q), and £1 is known 
as the dynamical phased varying from the value flLeft 
at the left shoulder of the diffraction peak to the value 
Qmght at the right shoulder; both values will be given 
latter on. 

To calculate the positioning of the standing waves dur- 
ing the rocking curve, the incident Dj and reflected Dr 
waves at a given depth h are approximated to 

D 7 = D (h)v e^M-feo-rO ( 9 ) 

and 

D R = R N {8)D {h)v e ^(u;*-fe-T-) ( 1Q ) 

where Vq and v are their oscillation directions. A com- 
parable reflection coefficient at all depths is assumed, i.e. 
Rn(9) do not depend on h. All dependence with depth 
is accounted for Do(h), which is a smooth function with 
very small variation over the lattice period. 

The standing wavefield is therefore D$w — Dj + Dr 
whose time-average intensity at any depth is proportional 
to 



Isw = ' = 1 + \Rn{0)\ 2 + 2\Rn(0)\v ■ «cos$, 

providing antinodes (maxima of intensity) at 

$ = 5 H + « - 2ttSH ■ r = 2mn (12) 
for every integer value of m. 
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For sake of simplicity, a single atomic layer per lattice 
plane period is assumed so that H = z/d, 

F H =e iS ^ d Y,fn , 

71 

and Sh = S2itzo / d where zq gives the atomic layer posi- 
tion in the crystal unit cell. If zq = 0, the unit cell origin 
would fall on top of the atomic layer. 
From Eq. (|12[1 we have the antinodes at 

z A (n) = z + Sd(n/27r~m) (13) 
as a function of the dynamical phase f2, and 



Az A = z A (tt Rlght ) - z A (Q Le f t ) = S (14) 

7T Z 

as the displacement of the antinodes during the crystal's 
rocking curve. Such displacement must be invariant re- 
garding the choice of signal S and, therefore, it can be 
assumed that 




FIG. 2: Reflection R and R, and transmission T and T, coef- 
ficients for electromagnetic plane-waves in very thin and uni- 
form planes of low-absorbing matter. Stokes relations for thin 
film optics— provide conservation of the counting rate \R+T\ 2 
plus \R + T\ 2 , at detectors Di and D2 respectively, regard- 
ing the total impinging intensity from the identical sources Si 
and S2 ■ The top and bottom insets illustrate the phase delay 
tp = —SkA'B = —SkB'A of the transmitted waves across the 
plane thickness. 



An = n 
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-St: 



(15) 



since the displacement is experimentally observed^ as 
been Az A = —d/2. 

It is widely known from dynamical diffraction theory^ 
that r^Le/t = 7T and ^lm g ht = 0. However, a clear state- 
ment should be made to emphasize that this phase vari- 
ation is relative to the global phase signal S = +1, as 
demonstrated below. 



D. Dynamical phase shift versus global phase signal 

Coherent scattering and transmission of electromag- 
netic waves by very thin and uniform planes of low- 
absorbing matter are describable by reflection, R and 
R, and transmission, T and T, coefficients as depicted in 
Fig. 2. The condition 



in order to fulfill Eq. OH). S± 90° and 5 ± 90° arc the 
amount by which the phases of the reflected waves can 
differ from the phases of the transmitted ones, and <p and 
Cp are the phase delays across the plane, whose thickness 
d is comparable to the wavelength A. 

To determine the phase delays consider first only the 
source S\ in Fig. 2. As shown in the top inset, if D(A) 
is the incident wave at point A, D(B) = TD(A) is the 
transmitted wave at point B passing through point A' . 
Then, 



\T\D e 



iS(ujt — ka-rB) 



TD e 



iS(ujt — ko-r A i) 



where fco ■ r A' — ■ r A since r A and r A i stand for posi- 
tions on the same wavefront AA' , providing 



T = \T\ 



-iS{r B - 



ITI 



-iSkA'B 



ITI 



Since the incidence angle 9 is the same for both sources 
A'B = B'A and, consequently, 



RT* + R*T + RT* + R*T = 



(16) 



is required by energy conservation, which stipulates 
phase relationships between the reflected and transmit- 
ted waves, as usually obtained for laser beam splitters^ 
or, equivalently, by the Stokes relations for time re- 
versibility of wave's propagation 11 ;. 

Without loosing generality, these coefficients can be 
written as 



2tt , . 
ip — ip ~ — S—dsmi 
A 



(18) 



The X-ray reflectivity of a crystal, as given in Eq. , 
is obtained by stacking lattice planes whose reflection and 
transmission coefficients are written in the same format 
used above, Eq. I|17|l . Stacking in geometrical progression 
of N — 2™ planes (n — 1, 2, ... ) is the fastest way to go 
from a single lattice plane to the desired thickness Nd. 
It is possible by means of the recursive equation 



R = ±i\R\ e i(s+ v\ T= \T\e lv , 
R = ±i\R\ e i{Uv \ and f = |T| 



(17) 



Rn 
Rn 



Tn/2Tn/2 
1 - Rn/2Rn/2 



N/2 
N/2 



(19a) 
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FIG. 3: Simulated X-ray reflectivity curves |.R(#)| 2 , for crystal thickness Nd = 2.5/xm (gray curve) and Nd = 21mm (black 
curve, flat-top maximum equal to 1). The shift of the dynamical phase Q, Eq. is shown on both sides of the rocking curves 
for the two possible choices of global phase signal in Eq. 1181 1. i.e. S = ±1. A8 — 9 — 9Bra gg , 2dsin6s ra99 = A = 1.54A (B ragg s 
Law), d = 3.14A, \Ri \ = \Ri \ — 3.0 x 10~ 4 (nearly the value for the silicon 111 reflection) and 5h = 0. For simulation purposes 
some absorption had to be considered in the thick crystal case, then a linear absorption coefficient fi — asin#B ragg /d ~ 1cm -1 
was used, a is defined in Eq. Il20bll . 



and 



III. DISCUSSION ON PHASE INVARIANTS 



T N 



1 — Rm/oR 



i JV 72 

T 2 

± N/2 



(19b) 



whose derivation is very similar to the Airy's formula of 
the Fabry-Perot interferometer, and where 



Ri 



.r e X\C\d 



V c sin t 



F H 



and 



(1 
(1 



\Ri? 
\Ri? 



a) 1 ' 2 
a) 1 ' 2 



(20a) 



(20b) 



C = Vq ■ v, H = H, r e = 2.818 x 10 _5 A (classical elec- 
tron radius), and a stands for photoabsorption probabil- 
ity on each individual lattice plane, as discussed in details 
and compared to other diffraction theories elsewhere^. 

Fig. 3 shows the behavior of the dynamical phase 
across the reflectivity curve as a function of the phase sig- 
nal S in Eq. 118|) . As expected, flLeft — w and Cljught = 
for the S = +1 choice, but flLeft = and O night — tt for 
the other choice, i.e. S = — 1, which is also in agreement 
with Eq. l(T5)l. 



Phase determination is a fundamental problem in X- 
ray crystallography. Reflection intensities only provide 
|.Fjr| as experimental data input for p c (r) in Eq. J7J). 
Then, what phase 5h should be assigned to the structure 
factor Fh = \Fh\ exp(i<5y) of each crystal reflection? 

For decades, great effort has been dedicated in develop- 
ing and improving methods to estimate reflection phases, 
but since these are relative choice-of-origin values, the 
estimable quantities are in fact the differences among re- 
flection phases, better known as phase invariants^. For 
instance, consider three reflections G, H and L whose 
diffraction vectors fulfill the condition G = H + L. 
Hence, in the structure factor ratio 



FrFl 
F g 



FrrWF; 



\F G 



the triple phase $ = 5h + <5l — <$g is invariant regarding 
the choice of origin in Fig. 1. 

Three-beam diffraction experiments are sen- 

sitive to 'J. In azimuthal scan mode, the diffracted in- 
tensity of reflection G is modulated by the excitation of 
the reflection H during the crystal rotation <p around G. 
This intensity modulation is approximately given by 



5 



I^) = \D G \ 2 + \D HL ^)\ 2 + 

+2\D a ■ D HL (4>)\ cos(* + fi), (21) 

which explicitly depends on the dynamical phase Q of the 
reflection H. Dq and Dhl{4*) stand for the diffracted 
wavefields from reflection G, kept constant during the 
0-scan, and from the detour reflection H + L. 

Although the signal of iff+ft, changes with S, the invari- 
ance of I(4>) is assured by the cosine, an even function. 
However, both angles must stand for the same choice 
of S in Eqs. (@J and (JTHJ. Otherwise, the experimen- 
tal determination of '5 by </>-scan data analysis would be 
compromised. 

Besides the signal of "J/, its modulus will also change 
with S when the /' and /" corrections 1 of the atomic 
form factor — the so-called anomalous scattering correc- 
tions — are taken into account. According to the struc- 
ture factor expression in Eq. 

\F H (S)\e^^ =^2(f + f' + if") n e iS2 * H - r % (22) 
n 

and then |F ff (+l)| ? \F H (-1)\ and 6 H {+1) S H (-1) 
when /" 7^ 0. ft is unacceptable since these inequalities 
would imply in different diffracted beam intensities for 
each possible choice of wave representation, S — +1 or 
5 = — 1. To avoid such artificial fact, the /" correction 
must bear its depence with S, so that 

f = h + .f + iSf". (23) 



IV. CONCLUSIONS 

Structure factor, dynamical phase shift, phase invari- 
ants, and complex atomic form factor are values suscep- 
tible to the choice of global phase signal, as explicitly 
demonstrated here. When a measurable quantity de- 
pends on more than one of these values, they must be 
calculated for the same choice of phase signal. 

Recursive equations based on Stokes relations are ap- 
plied for calculating the reflectivity of low-absorbing crys- 
tals with a finite number of lattice planes. These re- 
cursive equations are extremely simple and, yet, capable 
of describing important features of the diffraction phe- 
nomenon such as kinematical diffraction (thin crystals), 
primary extinction (maximum reflectivity equal to 1), in- 
trinsic width of Bragg reflections in either kinematical or 
dynamical diffraction regimes, and the phase shift of the 
diffracted waves across the reflection domain. 
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